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Decay constants of D and B mesons are estimated within the framework 



> 

in 
cn 

[-^ ' of a heavy-quark approach using measured isospin mass splittings in the 

O ■ D, D*, and B states to isolate the electromagnetic hyperfine interaction 

: between quarks. The values fo = (262 ± 29) MeV and fs = (160 ± 17) 
MeV are obtained. Only experimental errors are given; possible theoreti- 

CL|| cal ambiguities, and suggestions for reducing them, are noted. 

\ The decay constants and fs of mesons containing a single heavy quark are of 

fundamental importance for the understanding of the strong interactions, since they 
^ ■ describe the behavior of a single light quark bound to a nearly static source of color. 

The constant fs is crucial for interpreting data on particle-antiparticle mixing in the 
neutral B meson system, and both constants are essential if one is to anticipate and 
interpret new signatures for CP violation. In this Letter we describe a method for 
determination of these constants which relies on the isospin splittings of the D, D*, 
and B mesons, and indicate what additional information will be necessary to reduce 
systematic errors to an acceptable level. 

The decay constant /m for a meson M is specified by the matrix element of the 
axial current between the one-particle state and the vacuum: (0|y4^|M(g)) = iq^fu- 
In a non-relativistic quark model it is related to |\1'(0)|, the wave function at the 
origin, by [1] 

/m = (12/Mm)^/'|^(0)| . (1) 

While we recognize the limitations of the nonrelativistic model and the relation (1), 
especially for D mesons [2-4], we seek independent information on |^E'(0)|. We find 
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it by comparing isospin splittings in pseudoscalar and vector meson multiplets [5,6]. 
These have recently been measured very precisely, with the results [7] 

5m(D) = M(D+) -M(D°) = (4.80 ±0.10 ±0.06) MeV ; (2) 

5m(D*) = M(D*+) - M(D*°) = (3.32 ± 0.08 ± 0.05) MeV ; (3) 
5m(D) - 5m(D*) = (1.48 ± 0.09 ± 0.05) MeV . (4) 
There are previously existing measurements of the mass difference 



Sm{B) = M{B°) - M{B- 



(2.0 ± 1.1 ± 0.3) MeV CLEO 85 [8] 

(-0.4 ± 0.6 ± 0.5) MeV CLEO 87 [9] 

(-0.9 ± 1.2 ± 0.5) MeV ARGUS [10] 

(0.12 ±0.58) MeV (Average) 



(5) 



Within certain assumptions, the values (2) - (5) allow one to determine 1^^(0)1 for 
the D and B meson systems, and hence to calculate fo and fs- 

Three sources of isospin mass sphttings in hadrons can be identified [11]. 

(a) The d-u mass difference leads to splittings both directly and through the QCD 
hyperfine interaction 

const, {aj • aj) 

^-C/QCD hfs = ■ (Dj 

rriimj 

Here (a^ ■ aj) = (—3, 1) for a total spin Sij = (0, 1) of the ij quark pair. The 
constant depends on the color representation of the quark pair and thus is different 
for mesons and baryons. One can estimate its value by comparing masses of states 
with different spins, such as nucleon and A, or D and D*. Fits to hadron masses 
based on constituent quark masses, with account taken of the term (6), are remarkably 
successful [12,13]. In a nonrelativistic model with single-gluon exchange, the constant 
in (6) would be proportional to q;s|^'(0)P, where is the strong coupling constant. 

(b) Coulomb interactions among quarks lead to an energy shift AE'coui — QiQj{^/f)ij 
for each interaction between quarks of charges Qi and Qj. 

(c) The electromagnetic hyperfine interaction leads to an energy shift associated 
with the interaction of quarks i and j: 

. J, 27ra|^(0)pg,g,(a,-a,) 

^-C/e.m. hfs = 7, ■ i ' j 

iniimj 

In contrast to the strong hyperfine interaction, where the magnitude of the constant 
in (6) is uncertain a priori (though it can be measured experimentally), the term (7) 
holds the promise of providing information on |^(0)p. 

One can describe the isospin mass splittings in a meson in terms of three unknown 
parameters: m^- = x, the expectation value of 1/r, and the value of ^(0). Quark 
masses are assumed known from fits to meson mass spectra based on the constituent 
quark model [13]: mav = (m„ + md)/2 = 310 MeV, rUc = 1662 MeV, = 5 GeV. 
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The isospin mass splittings for charmed mesons are predicted to be [14]: 

where AMjj = 141.4 MeV is the D* — D mass difference (we use the average for the 
charged and neutral states). Here we have expanded the effect of 7^ m„ in the 
QCD hyperfine interaction term (6) to first order in x = — m„. For B mesons the 
corresponding expressions are 

where the B — B* mass splitting term AMb is about 46 MeV [15]. 

In the spirit of heavy-quark symmetry [16], we shall assume that the wave function 
of a light quark bound to a heavy quark is independent of both the identity and the 
spin of the heavy quark. With this ansatz, Eqs. (8) - (10) express the three measured 
quantities (2), (3), and (5) in terms of the parameters x = rrid — rriu, 

(1/r) ^ {l/r)D = {l/r)D* = (l/r)^ = (l/r)^* , (12) 

and 

|*(0)| = |*(0)|ij = |*(0)|d* = |*(0)|b = |*(0)|b* . (13) 
Solving, we find [17] 

a; = (1.29 ± 0.35) MeV ; a(l/r) = (3.60 ± 0.53) MeV ; 

|*(0)|2 = (11.3 ± 2.4) X IQ-^ GeV^ . (14) 

Taking Mm = [M(pseudoscalar) -|-3M(vector)]/4 in Eq. (1), we obtain the values (cf. 
U = 132 MeV, fK = 161 MeV) 

fo = (262 ± 28) MeV ; fe = (160 ± 17) MeV . (15) 

The value of fn is compatible with the experimental upper hmit [18] fo < 290 MeV 
(90% C.I.). It is in the range of recent theoretical estimates (see, e.g., [3,5]). The value 
of fs is compatible with analyses oi B — B mixing in the context of information about 

elements of the Cabibbo-Kobayashi-Maskawa matrix [19], but lies below predictions 
of lattice QCD [3]. These, as well as the work of Ref. [4], suggest that there are 
important corrections of order l/mg, where Q is a heavy quark, to the relation (1). 
The parameters (14) imply that 

8m{B*) = (0.08 ± 0.38) MeV . (16) 
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A measurement of dni{B*) with an experimental error comparable to that in (16) 
would be desirable as a check of our assumptions. 

The errors in Eq. (15) are purely experimental. There are both QCD and 0(l/mQ) 
corrections to the nonrelativistic relations (1) and (7). An example of the first is that 



the ratio fs/fD = y M^/Mb implied by the use of identical B and D wave functions 
in (1) should be multiplied by a QCD correction of 1.11 [20], while l/mq corrections 
to (1) have been discussed in Ref. [4]. However, until a corresponding discussion for 
the electromagnetic hyperfine term (7) has been given [21], it does not make sense to 
incorporate only partial information on such corrections. Corrections to the ratio of 
the two terms are what we need, and probably make more sense. 

The errors on light quark masses used in extracting |^(0)p from the electromag- 
netic hyperfine interaction energy are probably about 20%, based on the spread in 
values obtained in various constituent-quark fits to mesons and baryons [13]. These 
errors lead to uncertainties in fu and /b comparable to those due to uncertainties in 
|*(0)p in Eq. (14). 

Another way of finding x = nid — rriu is to use fits to the isospin mass splittings in 
the nucleon, S, and S states, based on the effects (a) - (c) mentioned above [14]. An 
independent estimate of x allows us to obtain by comparing (4) with the difference 
between (8) and (9) (the a{l/r) contributions cancel in the heavy-quark limit) without 
assuming that wave functions are identical in D and B systems. However, we shall 
see that a much larger value of x arises than one finds in Eq. (14). This provides 
us with an idea of systematic errors inherent in the use of the nonrelativistic quark 
model, and emphasizes the importance of knowing the full set of corrections to the 
ratio of (1) and (7). 

The QCD hyperfine interaction in baryons may be obtained from the mass dif- 
ference Abaryon = M(A) - M{N) = (1238 - 938) = 300 MeV. Neglecting fiavor- 
symmetry-breaking effects in wave functions but keeping them in quark masses, we 
define 



,1, 27ra,^,^, ,2 mo 



where mg = 363 MeV is the average u and d quark mass found in a fit to baryon 
masses, while m,, = 538 MeV in the same fit [13]. The expectation values of the ai-aj 
terms in (6) and (7) may be evaluated using quark model wave functions [13], while 
pairwise Coulomb interactions of quarks may be summed to obtain the coefficients of 
the term y. To first order in isospin-splitting terms, one can then write 

M(n) - M(p) =(l- ^^^) x--y + -z = 1.293 MeV ; (18) 
V 3mo / 3 3 



M(E-) - M(E+) = f 2 - + f^^^) x+ly+(^r + 

\ 3mo 3ms / 3 V3 

- (8.07 ±0.09) MeV ; (19) 
M(E+) -t-M(S-) -2M(E°) =1/-;^= (1.70 ±0.12) MeV ; (20) 
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1 + — ) x + + -r^ = (6.4 ± 0.6) MeV . (21) 
3ms / 3 3 

We have used baryon masses from Ref. [22]. Combining (18) and (20), we find (1 — 
[Abaiyon/3mo])a; = (1.86 ± 0.04) McV, or a; = (2.57 ± 0.06) MeV. The corresponding 
central values of y and z are 3.06 and 1.36 MeV, respectively. Comparing (4) with 
the difference between (8) and (9), we then find |*(0)|^(,) = (3.9 d= 1.3) x 10"^ GeV^ 
and fo = (154 ± 26) MeV for m^y = 310 MeV in (8) and (9). The comparison of 
this result (based on x from baryon masses) with our previous ones (14) and (15) 
(based on assumptions motivated by heavy-quark symmetry) shows that these two 
independent ways of deriving values for x give very different answers for Jo- 
lt is possible that one should not use the same values of x in mesons and baryons. 
However, unless one is prepared to make the assumptions leading to (14) and (15), 
there is not enough independent information available on x Prom mesons alone. The 
K and K* systems involve wave functions sufficiently different from the D and D* 
that their isospin splittings do not provide a reliable estimate of x. 

We have shown that the values fo = (262 ± 29) MeV and fs = (160 ± 17) MeV 
can be obtained from isospin mass splittings in the D, D*, and B meson systems, if 
the nonrelativistic formulae (1) and (7) relating the wave function at the origin to the 
decay constant and to the hyperfine interaction are valid, and if the wave functions 
are identical in the three systems. 

Our use of heavy quark symmetry in the present context can be checked not only 
by measurement of fo (via detection of the decay D — > /x+z/) , but also by verification 
of the very small predicted isospin splitting in the B* system. 

One might expect the relation (1) to be better satisfied for B mesons than for D 
mesons. A reliable estimate of x and a measurement of the difference between isospin 
splittings in the pseudoscalar and vector B mesons would permit an independent 
determination of fs using the present method. 

Considerable systematic error comes from uncertainty in constituent-quark masses. 
Just the uncertainty in whether to use light-quark masses from fits to mesons or to 
baryons introduces an additional error of about 10% in the decay constants. A much 
more serious problem is that finding x in another model-dependent manner using fits 
to the isospin mass splittings in the baryons gives a value of x about twice as large 
as that required to explain the isospin sphttings in the D, D*, and B systems. This 
potential inconsistency might be resolved once the full set of QCD and 0{l/mQ) 
corrections to the ratio of (1) and (7) have been calculated. 

This work was supported in part by the U.S. Department of Energy under Grant 
No. DE FG02 90ER 40560, and by the National Science Foundation. 
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